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Abstract 

In this paper thermal conductivity and thermal diffusivity of a two layer 
system is examined from the theoretical point of view. We use the one di- 
mensional heat diffusion equation with the appropriate solution in each layer 
and boundary conditions at the interfaces to calculate the heat transport in 
this bounded system. We also consider the heat flux at the surface of the 
samle as boundary condition instead of using a fixed tempertaure. From this, 
we obtain an expression for the efective thermal diffusivity of the composite 
sample in terms of the thermal diffusivity of its constituent materials whithout 

any approximations. 
44.30,+v, 44.10.+i 
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I. INTRODUCTION 



During the last decade, several methods have been developed to determine thermal dif- 
fusivities and conductivities with high precision by means of photothermal effects.0 The 
most widely used method is based on the photoacoustic effect. The principle of the effect 
is that when a sample in a closed cell is illuminated by light modulated or chopped at au- 
dio frquencies, an acoustic signal is produced. The application of the photoacoustic effects 
to the measurement of thermal diffusivities for thin films has been made for Adams and 
Kirkbright.! They have used the method called rear-surface illumination. When the rear 
surface of a sample is illuminated with the chopped light beam, heat oscillations gener- 
ated therein propagate from the surface into the sample. Pressure oscillations of the same 
frquency are induced in the gas chamber by the temperature oscillations at the surface in- 
terface between the sample and the gas, where they can be detected by a microphone. The 
photoacoustic signals obtained have certain phase shift relative to the signal detected from 
the front-surface excitation. Besides that, since the phase shift of the signal does not depend 
on the optical properties of the sample, it is simpler to extract information on the thermal 
diffusivity from the experimental results. Phase shift lag masurements are then most suit- 
able for determining the thermal diffusivity of the material. Charpentier et al.S made an 
analysis of the pressure variations considering rear and front-surface excitation and gave a 
formula relating the lag shift of the photoacoustical signal using high modulation frequency. 
The refinement of the theory, in terms of the Rosencwaig and Gersho theoryi was given by 
Pessoa et al.i , they showed the relative phase lag exhibits no explicit dependence on the 
absorbed power and surface conditions so that a single modulation frequency measurement 
is sufficient to determine the thermal diffusivity. 

In recent years, has been some interest in study the thermal characterization of two-layer 
system of variable thickness using the photoacoustic effect. Tominaga and ItoEl used the 
Rosencwaig- Gersho model to a two-layer system under rear illumination and looking at the 
phase angle behaviour as a function of the modulation frequency. They showed that, a high 
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modulation frequencies, the rear-illumination phase angle depends upon a critical frequency 
above which one of the materials becomes thermally thick. From the analogy between 
thermal and electrical resistances widely used in heat-transfer problems,! Mansanares et 
al.@ calculated the effective thermal diffusivity of the two layer system as a function of the 
filling fraction of the composite sample and the ratio of the thermal conductivities of each 
material. Recently Christofides and Seas! made an extension of the theoretical model on 
photopyroelectric spectroscopy of solidsB to investigate the optical and thermal properties 
of a two-layer sample. They examined the case where the optical absorption coefficients of 
the substrate and film vary in different ways, several computer simulations were performed 
in order to examine the correctness of this model for a wide range of wavelengths and 
modulation frequencies. 

In the present paper, we take a different approach than previous investigations on the 
heat transport in bounded systems. We restrict our analysis to the case of the solution of 
heat transport equation only considering the continuity of the heat flux at the interface of 
the two-layer system. In addition, we only take into account the decreasing exponential term 
of the dynamical part of the temperature fluctuation in each layer which gives the physical 
solution of the heat flux. It is shown, that if the ratio of the square root of the diffusivity is 
equal to the ratio of the thermal conductivity of each layer, then the temperature distribution 
is continuous at the interface, otherwise the temperature is discontinuous. In both cases, 
the effective diffusivity of the two-layer system is the same. 



II. ONE-LAYER SYSTEM 

It is well known, that heat transport in solids is carried, by various quasiparticles (elec- 
trons, holes, phonons, magnons, plasmons, etc). Frequently the interactions between these 
quasiparticles are such that each of these subsystems can have its own temperature and the 
physical conditions at the boundary of the sample can be formulated separately for each 
temperature. For example, the physical conditions resulting in heat transport by electrons 
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and phonons are given by Granovski and Gurevicb0 and those for the transport by electrons 



and magnons are given in Ref. |T^, where the appropriate boundary conditions are formu- 
lated. The theory of heat conduction in solids was developed by Gurevich and Kaganov0 
using the two-temperature approximation. They showed that, in general, the temperature 
of carriers and phonons in anisotropic semiconductors are unequal even in the interior of a 
bulk sample. The expression for the temperature distribution for electrons and phonons was 
obtained in bounded semiconductors using the adiabatic boundary conditions. It is showed 
that the accepted assumptions about the constancy of the temperature gradient is only valid 
under certain limits. 14 

We restrict ourselves for defmiteness to the case when the quasiparticle systems are 
electrons and phonons. Let T p the characteristic phonon temperature. Then the momentum 
q oc — , s is the sound velocity, sets a limit to the phase space volume occupied by phonons. 
Hereinafter, we shall consider two limiting cases:E3 

(i) Short wavelength (SW) phonons occupy a large volume in phase space. This case is 
represented by the the inequality 

2 P « (1) 



where p is the average electron momentum, namely p oc ^2mT e for nondegenerate electron 
gas, and p — Pf , the Fermi momentum, for degenerate electron gas, here T e is the electron 
temperature. 

(ii) Large wavelength (LW) phonons occupy a small volume in the phase space. This is 
the case when electrons interact with all phonons, 

2p > -A (2) 

From Boltzman equation, it is known that the degree of nonequilibrium phonons is 
determinated by the relationship between phonon-electron (v pe ) and phonon-phonon (v pp ) 
relaxation frequency, v pe determines the degree of phonon disturbance by the electrons and 
it decreases rapidly tending to zero for q > 2p, while the other, u pp , describes the tendency 



of phonons to come to equilibrium as a result of energy distribution. Therefore, in the region 
q > 2p, only the following inequality holds 

v w > Vpe, q> 2p. (3) 

In the region q < 2p, both the inequality 

v vv > v pe , q <2p (4) 

can hold as well as the reverse inequality 

v w < Vpe, Q> 2p. (5) 

In the limit v pp ^> i/ pe , the phonon-phonon interations are more frequent than phonon- 
electron collisions and more efficient in terms of energy relaxation than energy transfer from 
the electron to the phonon subsystem. It should be noted that the electron and phonon 
subsystems, generally speaking, cannot be characterized by a single temperature. Therefore, 
steady state heat conduction can be described by the following system of equations 

divQ e = P ep (T e - T p ), divQp = -P pe {T e - T p ) (6) 

the term P ep {T e — T p ) describes the transfer of heat between electrons and phonons. Here 
P ep is a parameter proportional to the frequency of electron-phonon collisions (P ep = P pe ) 
and the heat flux of electron Q e and phonon Q p subsystems are described by the usuall 
relationships 

Q e = -£; e divT e , Q p = -£; p divTp (7) 

where k e (k p ) is the electron (phonon) thermal conductivity. If in addition, we consider 
large specimens such that the dimension I 3> L*, where L* is the scale of the electron- 
phonon energy interaction, referred to as the cooling length then the temperature of the 
two subsystems are equal i.e. T e = T p = T.0 In this case we obtain, after sum Eqs.(6), the 
following equation 
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divQ = (8) 

where Q = Q e + Q p is the total heat flux carried by electrons and phonons. 

Now, consider the situation represented by the inequality of Eq.(5). This is the case 
when phonon-phonon collisions alone can not bring the phonon subsystem to an internal 
equilibrium. If Eq.(l) is true i.e. LW phonons occupy a much smaller phase volume than SW 
phonons then, the SW system has enough time to redistribute the energy recived from the 
LW phonons between its consituents quasiparticles. As a result, the distribution function of 
SW phonons becomes Planckian (see Ref. [15]). Then the electron- LW phonon interactions 
relax its energy more efficient that the phonon subsystem and the LW phonons emitted 
by electrons of temperature T e are characterized by the same temperature T e = T p . In 
this situation we also have two different subsystems, one corresponds to the SW phonons 
with temperature T p w and the other one corresponds to electron and LW phonon with a 
characteristic temperature T e = T p w and they satisfay the following heat transport equations 

div(Q e + Q£ w ) = -P pp (T e - T p LW ), divQ p sw = P pp {T e - T p LW ) (9) 

where the term P pp (T e — T p w ) represents the transfer of heat from LW phonons to SW 
phonons and P pp is calculated in Ref. [T^. If the size of the sample is greater than the cooling 
length (7 ^> Lq) of this system then we obtain the condition T e = T p w = T and after sum 
Eqs.(9), the total heat flux satisfies divQ = 0. 

When v pe 3> v pp and Eq. (2) holds, we have only LW phonons and they interact efficiently 
with electrons and in this case all the quasiparticles system have the same temperature 
T e = T p = T and the equation for heat flux is similar to Eq.(8). 

From this brief discussion, we have shown that under certain conditions on the realaxation 
frequency of electron and phonon subsystems and size of the sample the total system can 
be described by the same temperature T and for the total heat flux divQ = 0. 

So far, we have only described the static contribution of the heat transport i. e. the 
heat flux is independent of time. However, in the photoacoustic experiments, the incident 
radiation is modulated on time by the chopper and in this case it is necessary to consider the 
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dynamic contribution on the heat transport in the electron and phonon system. It is worth 
to mention that even the external perturbation depends on time, the dynamic contribution 
of the heat flux is only important when the frequency u of the incident electromagnetic field 
is of the same order of the characteristic relaxation energy frequency of the electronic system 
i/ e .ll3 If ijj ^> v e (high frequency limit), in this case the electron temperature cannot follows 
the variation of the field and assumes an average value. Since phonons can recive energy 
only from electrons, the phonon temperature should also remain constant as function of 
time. Otherwise, if uo <C u e , (low frequency limit) the variation of the electron and phonon 
temperature is quasi static. That means, the static quasiparticle temperature oscillates with 
the same frquency of the radiation. 

In Eq.(8), we are not taking into account the distributed heat source resulting from the 
ligth absorption. 

Assuming that the temperature of the electron and phonon systems are equal, which 
is usually the condition in most of the photothermal experiments the equation of heat 
conduction in solids which is valid for u < v e can be written as 

^|^ = «V 2 T(r,t). (10) 

In this equation, we are considering that the variation of the temperature as function of 
x is such that the heat conductivity k is independent of the coordinates, otherwise, we have 
to solve a non-liner heat equation. Here the diffusivity a is given as a = where p is the 
density and c is the specific heat of the sample. Consider the photoacoustic cell geometry 
for the heat transmission configuration shown schematically in Fig. la. The temperature 
fluctuation is obtained from the solution of Eq.(10) in one dimension. The solution T(x,t) 
should be supplemented by boundary conditions at x — 0. In the photoacoustic experiments, 
the most common mechanism to produce thermal waves is the absorptin by the sample of 
an intensity moduladed light beam with frequency modulation uj < v £ . It is clear that fixed 
the intensity of the radiation, the light-into-heat conversion at the surface of the sample can 
be written as 
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Q{x,t)\ 



x=0 



Q + AQe 



(11) 



where Q is proportional to the intensity of high frequency light (Q ^> v £ ) and the other 
term represents the modulation of this light. The temperature is not used as boundary 
condition because it is usually an unknown parameter in the experiments and besides that it 
is necessary to know the temperature on both surfaces. It is only important in thermoelectric 
phenomenon in semiconductors when the specification of the temperatures on the surfaces 
of the sample must to be known. 

The general solution of the heat diffusion equation for one-layer system is given by 



The parameter a is determined by forcing Eq.(12) to satisfies Eq.(lO) for one dimensional 



heat flux and is equal to a — (1 + i)yu/2a and using the boundary condition at x — 0, the 
constants 7\ and T 2 are given by 



In arriving to Eq.(12), we assume that the sample is optically opaque to the incident 
light (i.e. all the incident light is absorbed at the surface), here T is a constant which 
cannot be determined from this boundary conditions and it is not important in obtain the 
physical results. It is worth to mention that the increasing exponential term with distance 
in the dynamical part of Eq.(12) has not been considered because this term represents a 
macroscopic heat flux from the lower to higher temperature region (heat flux cannot be 
reflected) . 

Once we know the temperature distribution in the sample, we can assume the acoustic 
piston model for evaluating the pressure fluctuation in the cell. Acording to this model,§ the 
ocillating component of the temperature attenuates rapidly to zero with increases distance 
from the sample surface. The thin gas boundary layer at the interface is then to be acting 
as a vibrating piston. The displacement of this piston is estimated using the ideal gas law 



T(x,t) = T + T 1 x + T 2 e 



(12) 





(13) 
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(PV = NksT) for the boundary layer. As a result of this gas piston oscillation, a pressure 
fluctuation 5P(t) is produced in the cell and is given by 

«w = tM (£) " 2 K< - ^ + - i/P - (£) 1/2 <"> 

where P an d are the ambient pressure and temperature respectively. Then Eq.(14) may 
be evaluated for the magnitude and phase of the acoustic pressure wave produced in the cell 
by the photoacoustic effect. 



III. TWO-LAYER SYSTEM 

Let us consider the two-layer system shown schematically in Fig. (lb) consisting of a 
material 1 of tickness l\ and of material 2 of tickness l 2 , both having the same cross section. 
Let L — li + l 2 denote the total sample thickness, the thermal diffusivity and ki the 
thermal conductivity of the material i (i — 1, 2). The system of heat diffusion equations 
describing the heat transfer through the various layers of the one dimensional photoacoustic 
configuration are given by 

9Ti d 2 Ti 

-m= a ^- (15) 

The boundary conditions of the thermal diffusion equation (15) are obtained from the 
requirement of the heat flux continuity at the interfaces of the two materials and Eq.(ll). 
The solutions of the heat transport equations can be written as 

Txfot) = T - 9- x +^Q-(l-i) (^V /2 e-*-i*, < x < h, (16) 



ki 2ki V uj 

T 2 (x,t)=6o + 9i(x-l 1 ) + 9 2 exp[iu}t-cr2(x-li)}, h < x < L, (17) 



where 



9i = (18) 



2 
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1 - z) exp 



2k2 \ uj 

with 



(19) 



« = (£lT<- 1 + < >- < 20 > 

Comparing temperature distribution Eq.(17) with Eq.(12) at x = L, we can write the 
effective thermal diffusivity and conductivity of the two-layer system as 

_L = J^ + A_ (21) 

and 



k ko 



(22) 



or 



012 ' h +-%=). (23) 



k k<2 \yfOL\ y/02 

It is worth to mention that effective thermal parameters Eqs.(21) to (23) have been 
obtained whithout any approximations about the thermal thin and thickness materials]! 
analogy between thermal and electrical resistances used in heat transfer problems! or some 
critical frequencies above which one of the layers becomes thermally thick! and continuity 
of the temperature distribution at any interface. 

However, if the thermal parameters of each layer satisfy 



«i V« 2 (24) 



k\ &2 

and using Eq.(21),the effective thermal conductivity can be written as 

H + l (25) 

Then, from Eqs.(16) and (17) we obtain that Ti(li,t) = T 2 (li,t) i.e. the temperature is 
continuous at the interface x = l\, otherwise, in general, the temperature will be discontinous 
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at the interface of the two materials and effective thermal conductivity and diffiusivity 
are given by Eqs.(21) and (23). This dicontinuity in the temperature can be expressed 
mathematicaly as follows: 

Q(M)| a=ti = gm-fc g^ +£) ~ Tl(/l ~ £) = V (T 2 - T x ) (26) 

where r] is the surface thermal conductivity at the interface.0 Note that when r] goes 
to infinite, since Q is finite, the temperature distribution must continuous at the interface 
and for the surface thermal conductivity finite, in general, we have that Ti(x,t)\ x=li ^ 
T2(x,t)\ x=h . 



IV. CONCLUSIONS 

A theoretical analysis of the photoacustic effect on a two layer sample has been studied. 
Using the appropriate boundary conditions, we obtain the effective thermal diffusivity and 
thermal conductivity of the two-layer system whithout any approximation on the thermal 
parameters. The continuity or discontinuity of the temperature at the interface of the 
two-layer system depends on the relationship of the thermal parameters of both layers. In 
general, the heat flux is defined as the product of the thermal surface conductivity and the 
diference of the temperatures at the interface. 

Mansanares et.al.i demonstrated the usefulness of a single modulation frequency method 
for measuring the thermal diffusivity of a solid samples. The method consists of measure 
the relative phase between the rear-surface illumination and the front-surface illumination. 
Using the thermal diffusion model of Rosencwaig and Gersho§ for the production of the 
photoacustic signal, the ratio of the signal amplitud and the phase lag for rear and front 
surface illumination are given as function of the sample thickness and the sample thermal 
diffusion coefficient. The theory for the relative phase lag will be studied using our model in a 
forthcoming paper. Finally, its important to mention that our model is valid for modulated 
frequency of the incident ligth u, of the same order that the frequency of the relaxation 
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energy between the quasiparticle systems v e . In the limit, uo ^> v £ the system can not 
respond to this external perturbation, therefore, the average in time of the dynamical part 
of the heat flux is neglegible and the transferred heat is only static. For u <C u £ , the heat 
flux is quasi-static i.e. dT/dt = 0.0 Then the solution of Eq.(10) for one layer is given 
by T(x,t) = C\ + C2X and from Eq.(ll), we finally obtain the solution in this regime as 
T(x,t) = (C[ + C^xje 1 ^, this represents an oscillation of the temperature distribution in 
every point of the layer. 
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FIGURES 

FIG. 1. Geometry for (a) one-layer system and (b) two-layer system. 
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